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ABSTRACT 

We discuss a simple model for the formation of a supernova spectral line by 
resonance scattering in the presence of external illumination of the line-forming 
region by light from circumstellar interaction (toplighting) . The simple model 
provides a clear understanding of the most conspicuous toplighting effect: a 
rescaling or, as we prefer, a "muting" of the line profile relative to the continuum. 
This effect would be present in more realistic models, but would be harder to 
isolate. An analytic expression for a muting factor for a P-Cygni line is derived 
that depends on the ratio E of the toplighting specific intensity to the specific 
intensity from the supernova photosphere. If £^ < 1, the line profile is reduced 
in scale or "muted". If £^ = 1, the hne profile vanishes altogether. If £^ > 1, 
the line profile flips vertically: then having an absorption component near the 
observer- frame line center wavelength and a blueshifted emission component. 

Subject headings: radiative transfer — supernovae: general — supernovae: 
individual: SN 1998S. 

1. Introduction 

In a simple but useful model of spectral line formation during the photospheric 
phases of a supernova, a line forms by resonance scattering in a homologously expanding 
atmosphere above a sharp photosphere. An unblended line has a P-Cygni profile with 
an emission feature near the observer-frame line center wavelength and a blueshifted 
absorption. Synthetic spectra calculated on the basis of this simple model have been 
found to fit the observed spectra of most supernovae reasonably well, and to be useful for 
establishing constraints on the composition structure of the ejected matter. For a detailed 
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discussion of the model, see Jeffery & Branch (1990), who also provided an atlas of line 
profiles, and compared synthetic spectra to numerous observed spectra of the Type II 
SN 1987A. Recent studies based on the simple model, using the SYNOW supernova 
synthetic-spectrum code, include Millard et al. (1999) on the Type Ic SN 19941 and Hatano 
et al. (1999) on the Type la SN 1994D. 

In the bright Type Iln SN 1998S we have encountered an event for which the above 
simple model fails. The supernova SN 1998S has been observed extensively with the Hubble 
Space Telescope (HST) as a target of opportunity by the Supernova INtensive Study 
(SINS) group (Garnavich et al. 1999, Lentz et al. 1999, and Blaylock et al. 1999, all in 
preparation), as well as from the ground (Leonard et al. 1999). The observed spectra, 
especially at short wavelengths and early epochs, contain numerous narrow (< 500 km 
s~^) absorption and emission features that formed in circumstellar matter. The optical 
and ultraviolet spectra also contain broad (~ 5000-10,000 km s~^) features that formed in 
the supernova ejecta. All of the broad features show little contrast with the continuum, 
especially at early times. It is not plausible that the radial dependence of the optical depths 
of all of the supernova lines should be such that the line profiles come out to be shallow, 
and in synthetic spectrum calculations with the SYNOW code in its simplest form we 
cannot account for the relative strengths of the spectral features of SN 1998S. The presence 
of the circumstellar features suggests that the broad supernova features were being affected 
by light from the region of circumstellar interaction: i.e., the line formation region was 
illuminated not only from below by light from the photosphere, but also from above by light 
from the circumstellar interaction. We refer to this external illumination of the supernova 
line-forming region as "toplighting" . 

In this paper we present a simple model of resonance-scattering line formation in the 
presence of toplighting. This simple model provides a clear understanding of the most 
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conspicuous toplighting effect: a rescaling or, as we prefer to call it, a "muting" of the line 
profiles relative to the continuum. This effect would be present in more realistic models, 
but would be harder to isolate. 

Section 2 presents the model. The emergent specific intensities are derived in § 3. The 
line profile and a muting factor are derived in § 4. In § 5, we offer a picturesque description 
of radiative transfer in the model atmosphere. A final discussion appears in § 6. 

2. The Model 

Suppose that supernova line formation takes place by isotropic resonance scattering 
in a spherically symmetric atmosphere that has a sharp photosphere with radius i?ph as 
an inner boundary and a circumstellar-interaction region (CSIR) as an outer boundary. 
Assume that the CSIR is a shell of zero physical width and optical depth, and let its 
radius be i?cs (see Fig. 1): R^s > Rph in all cases, of course. For simphcity, assume that 
the photosphere emits an outward angle-independent specific intensity /ph, that the CSIR 
emits isotropically with specific intensity Ics, and that /ph and Ics are constant over the 
wavelength interval of interest. When 7cs is set to zero, we have what we will call the 
standard case. Nonzero Ics gives the toplighting case. 

We assume that the atmosphere is in homologous expansion. In homologous expansion 
the radius of a matter element r is given by 

r^vt , (1) 

where v is the constant radial velocity of the matter element and t is the time since 
explosion which is assumed sufficiently large that initial radii are negligible. 

Line formation is treated by the (non-relativistic) Sobolev method (e.g., Rybicki 
& Hummer 1978; Jeffery & Branch 1990; Jeffery & Vaughan 1999) in which the line profile 
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in the emergent flux spectrum depends on the radial behavior of the Sobolev line optical 
depth T and the line source function S (which is a pure resonance scattering source function 
following our earlier assumption) . We consider only isolated (unblended) line formation and 
do not include any continuous opacity in the atmosphere. 

We note for homologous expansion that the resonance surfaces for observer-directed 
beams are planes perpendicular to the line-of-sight. (A resonance surface is the locus of 
points on which beams are Doppler shifted into resonance with a line in an atmosphere 
with a continuously varying velocity field.) If wc take z to be the line-of-sight coordinate 
with the positive direction toward the observer, the resonance plane for a wavelength shift 
AA from the line center wavelength Aq is at 

z^-—ct, (2) 

where z/t is the plane's velocity in the ^-direction. Blueshifts give positive z planes and 
redshifts, negative z planes. 



3. The Emergent Specific Intensities 

For resonance scattering the source function of a line is just equal to the mean intensity. 
In our model the line source function is 

S{r) = WI^^ + [l-W]I^, (3) 

where W is the usual geometrical dilution factor, 




(e.g., Mihalas 1978, p. 120). The first term in equation (3) accounts for radiation from the 
photosphere and the second, for radiation from the CSIR. 
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Consider a line-of-sight specific intensity beam that does not intersect the photosphere 
(i.e., one that has impact parameter p > Rp^ in the standard p, z coordinate system) and 
that has a wavelength in the observer's frame that differs from the line center wavelength 
by A A. In a tophghting case such a beam originates from the CSIR as shown in Figure 1. 
Prom the Sobolev method, the emergent specific intensity of such a beam can be seen to be 

lAxip > ^Ph) = ^cs e-" + [ly/ph + (1 - W%s] (1 - e-^) + 7cs , (5) 

where, of course, W and r must be evaluated at the Sobolev resonance point for the AA 
value under consideration. Similarly, the emergent specific intensity of a line-of-sight beam 
that originates from the photosphere (i.e., has p < Rp^) can be seen to be 

lAxip < Rph) = /ph e-^ + [WIp^ + (1 - W)/es] (1 - e-^ + /es ■ (6) 

Equations (5) and (6) can be rearranged into the convenient expressions 

Iax{p > i^ph) = (/ph - /cs)W^(l - e-^) + 2/es (7) 

and 

Iax(p < i?ph) = (/ph - /cs)e-" + (/ph - Ics)W(l - e-^) + 27,, . (8) 

If /cs = 0, these expressions reduce to the standard or non-toplighting case Sobolev 
expressions for emergent specific intensity. If the resonance point is not in the line-forming 
region (i.e., it occurs at r < /?ph or r > /?cs, or it is in the occulted region from which 
no beam can reach the observer), then r is just set to zero in the expressions which then 
reduce to the continuum expressions 

/aa(p > /?ph) = 2/,, (9) 

and 

/AA(p</?ph)=/ph + /cs- (10) 
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The emission component of a P-Cygni line is largely due to /aa(p > Rph) beams and 
the absorption component, to Iax{p < Rph) beams. To see how going from the standard 
to the toplighting case affects the components it is best to consider equations (7) and (6). 
Prom equation (7), we see that toplighting tends to reduce the line emission by changing 
/ph to /ph — Ics- ^Prom equation (6), we see that toplighting tends to fill in the absorption 
trough (caused by the e~'^ factor of the /phe^"^ term) by adding a positive term to the source 
function. Both emission component and absorption trough are reduced relative to the 
continuum by the addition of continuum terms: 27cs in the first case and Ics in the second. 
Thus we can see (at least for /ph > hs) that adding toplighting to an atmosphere is hkely 
to reduce the relative size of line components (i.e., to mute them). When /ph = Ics, the line 
components should vanish altogether as comparing equations (7)-(8) and equations (9)-(10) 
shows. In § 4, we give a definite analytical analysis of the effect of toplighting on line profile 
formation and confirm the muting effect. 

4. The Line Profile and the Muting Factor 

The flux proflle seen by a distant observer is obtained from 

rRcs 

Fax = 27r / dpplAxip) , (11) 

where the integration is over impact parameter. (Actually, the quantity in equation (11) 
divided by the square of the distance to the observer is the flux measured by the observer. 
But for brevity here and below we just call this quantity flux.) We can obtain semi-analytic 
expressions for the flux in the standard and toplighting cases by breaking the integration 
for the flux into components. These semi-analytic expressions then allow a fully analytic 
formula for a muting factor which describes the muting effect of toplighting. 

By inspection of equations (7)-(ll) with Ics = 0, we obtain the standard-case results 
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for flux in the continuum and in the hne. The continuum flux is 

F(con) = JphFo , (12) 

where 

>2 



The line flux is 



where 



Fo = nR' . (13) 



i^AA = /ph(Fl+F2) , (14) 



/ dppe-^ (15) 
Jo 



and 



F2 = 27r/ dppiy(l-e-^) . (16) 

Note that the Fi and F2 factors are dependent on z and therefore on wavelength. Also recall 
that T = for a resonance point outside of the atmosphere or in the occulted region. The 
contrast factor (i.e., relative difference of line flux from continuum flux) for the standard 
case is 

Faa - F(con) ^ Fi + F2 - Fp 

F(con) Fo ' ^ ^ 

Again from equations (7)-(ll), but with nonzero Ics, we obtain by inspection the 
toplighting-case results for flux in the continuum and line. The continuum flux is 

F*°P(con) = JphFo + 7es (2^0 - i^o) = (V - Ics) Fo + 2I^Go , (18) 

where 

Go = nRl . (19) 

The line flux is 

^AA = (^Ph - ^s) (i^i + F2) + 2I^Go . (20) 
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The 2IcsGq terms in equations (18) and (20) just account for the toplighting specific 
intensity beams aimed toward the observer from both the near and far hemisphere of the 
CSIR. These beams of course can interact with the hne and the photosphere and this 
interaction is accounted for in the other terms. Note that Go > -^o and Co > -Fi, where 
the equahties hold only in the degenerate case where i?cs = -Rph- Also note that Gq > F2, 
where the equality holds only in the degenerate case where both coefficients are zero: i.e., 
Res — Rph — 0. From these inequalities, it is clear that F*°P(con) and F^^ can never be less 
than zero: a physically obvious result, of course. 



The contrast factor in the toplighting case is 

^top _ irtop(con) ^ (/ph - /cs) {Fi + F2- Fq) ^ (1 - E) (Fi + F2 - Fp) 
Ftop(con) (7ph - /cs) Fp + 2I^Go (1 - E)Fo + 2EGp 



(21) 



where 



E=^ . (22) 



We define the muting factor m to be the ratio of the tophghting-case contrast factor to 
the standard-case one: 

1- E 

m = ^ . 23 

l-E + 2E{RjRpi,f 

Since all Fi and F2 factors have cancelled out, m is fully analytic and wavelength 
independent. The muting factor m is a monotonically decreasing function of E with a 
physical maximum of 1 at F = and with only one stationary point, a minimum at F = 00. 
The muting factor in fact goes to zero a,t E — 1 and becomes negative ior E > 1. Thus 
for F = 1 the P-Cygni profile of a line vanishes and for F > 1 the line flips: there is an 
absorption feature at the line center wavelength and a blueshifted emission feature. (Note 
that this flipped P Cygni proflle is not the same as the "inverse" P Cygni proflle that would 
be produced by a contracting rather than an expanding line-forming region.) A picturesque 
way of understanding a flipped P- Cygni line is given in § 5. 
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The minimum value of m at = oo is given by 

m{E = oo) = — 2 ■ (24) 

—1 + 2 {Res/ Rph) 

Note that m{E — oo) > —1 with the equahty holding only for Res /Rph — 1 which is the 
lower limit on Rcs/Rph- Thus 

|m| < 1 . (25) 

Since the absolute value of m is always less than or equal to 1, toplighting always mutes 
a line profile: hence our choice of "muting" for the name of the toplighting effect on line 
profiles. 

We can consider a couple simple examples of muting by toplighting. First, consider 
E — 1/2 and Res/ Rph — 2. Note that geometrically thin supernova atmospheres have not 
been identified, and thus Res/ Rph ^ 2 may be typical in real supernovae. With the given 
input values, m — 1/9 and the contrast factor of the line with respect to the continuum is 
reduced by this factor in going from a standard case to a toplighting case. The contrast 
factor of an absorption trough of a standard-case P-Cygni line has absolute lower limit of 
— 1 (i.e., zero flux). Thus, even if this lower-limit case were realized for a standard-case 
P-Cygni line, the toplighting counterpart absorption would have absorption trough depth 
of only 1/9 of the continuum level. 

Next consider E = oo and Res/ Rph = 2. Here m = —1/7. Since m is negative, the 
toplighting has produced a flipped P-Cygni profile. The standard-case P-Cygni absorption 
trough would be turned into a toplighting-case emission peak with an upper limit on the 
contrast factor of 1/7. 

Figure 2 shows examples of P-Cygni files with E values that effectively span much of 
the E parameter range. 

Since there is no upper limit on the contrast factor of a standard-case P-Cygni emission 
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peak, prima facie it seems that a negative muting factor with large absolute value could lead 
to negative observed flux in a corresponding toplighting-case absorption component. This 
does not happen of course. As we showed above, the observed flux in the toplighting case is 
never mathematically negative. For another point of view, consider the following argument. 
The standard-case P-Cygni line emission peak is largest for an opaque line (i.e., one with 
r — oo everywhere in the atmosphere) and large Rcs/Rph- For such a line, the emission 
peak contrast factor can only increase with increasing Rcs/Rph as ~ ln(i?cs/-Rph) (-Res 
considered here just as an outer boundary radius) (e.g., Jeffery & Branch 1990, p. 189), but 
the absolute value of the muting factor for large Rcs/Rph decreases for increasing Rcs/Rph 
like {Rcs/Rph)~^- Thus the muting factor always scales to prevent negative observed flux 
from arising mathematically. 

The ratio of the monochromatic luminosities of the CSIR and the photosphere (taken 
independently of each other) is 



where the factor of 2 accounts for the fact that both hemispheres of the CSIR region 
contribute to flux in any given direction. The ratio F, in fact, appears in the denominator 
of the muting factor formula, equation (23). If we use equation (26) to eliminate E from 
equation (23), we obtain 



Now m as a function of F monontonically decreases with F from m = latF = Otoa 
minimum m = — 1 + 2(i?cs/-Rph)^] at F = oo (the only stationary point). It is clear 
that F will have to be large in some sense in order to obtain strong muting. For the sake of 
deflniteness say m < 1/2 is "strong" muting. Then 




(26) 



2(i?cs/-Rph)^ — F 



(27) 



2{RjRphr-T + 2{RjRphn ' 




(28) 
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is required for strong muting. Since R^s/ Rph > 1 is required physically, a necessary, but not 
sufficient, condition for strong muting is F > 2/3. If, as suggested above, i?cs/-Rph ^ 2 for 
supernovae, then a necessary, but not sufficient, condition for strong muting in supernovae 
is r > 8/9. Consequently, only those supernovae whose monochromatic luminosities are 
strongly enhanced by circumstellar interaction will have line profiles that are strongly 
muted by toplighting. 



To complement the mathematical description of the radiative transfer in the standard 
and toplighting cases of our simple model, we present here a picturesque description. 

First consider the standard case. We use Res just as an outer boundary of the 
atmosphere in this case. If no line is present, the photons emitted by the photosphere just 
escape to infinity and the flux in any direction is just /ph-^o- this is the continuum emission. 
Adding a resonance-scattering line to the atmosphere has the overall effect of reducing the 
wavelength-integrated emission of the supernova in the wavelength interval that the line 
can affect. This interval is (Aq + AAmin, A -|- AAmax), where 



The reason for the loss in wavelength-integrated emission is that the line, scattering 
isotropically, will scatter some photons back to the photosphere where in our simple 
model they are simply absorbed. In a more realistic model, the photons absorbed by the 
photosphere are a feedback that helps determine the photospheric state. The particular 
observer-directed photons which are absorbed are those scattered toward the observer in 
the occulted region: they simply hit the photosphere as they head toward the observer. 

Because of the homologous expansion, photons continuously redshift in the comoving 



5. Picturesque Description 




and 




(29) 
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frame of the atmosphere. Thus, photons emitted by the photosphere at or redward of Aq 
escape the atmosphere without scattering. (Note formally they can scatter at Aq at the 
point of emission on the photosphere, but this effect is assumed accounted for in specifying 
/ph, the constant photospheric specific intensity.) Thus the line scatters photons that before 
scattering are blueward of Aq in observer-frame wavelength. The line scattering does not 
change a photon's comoving frame wavelength (not at all in our simple model and not 
significantly in reality) , but by changing its propagation direction relative to the matter flow 
it does change its observer-frame wavelength. Consider photons scattered toward a distant 
observer. If the scattering occurs from matter moving away from, mainly perpendicular to, 
or toward the observer, then the scattered photons in the observer frame have wavelengths 
redward of Aq, near Aq, or blueward of Aq, respectively. 

From the proceeding remarks, it is clear that the observer receives all the photons 
emitted by the photosphere into the line-of-sight at and redward of Aq plus the unscattered 
photons from blueward of Aq. In addition there are photons scattered into the line-of-sight 
by the line from the whole wavelength interval (Aq + AAmin, A + AAmax)- The scattering 
component is strongest near Aq where the resonance planes for scattering are largest and 
they come closest to (and even touch) the photosphere where the source function and 
usually the scattering opacity are largest. (Because of the nature of homologous expansion 
all the photons scattered toward the observer from a plane perpendicular to the line-of-sight 
have the same observer-frame wavelength: see § 2.) These planes are near z — 0. The 
scattering component grows progressively weaker as one moves away from near Aq: i.e., as 
the resonance planes get farther from z ^ 0. This scattering behavior results in the P-Cygni 
profile emission feature with its peak near Aq. 

The scattering of photons out of the line-of-sight from blueward of Aq causes a flux 
deficiency relative to the continuum: this is the P-Cygni profile blueshifted absorption. 
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As we argued above the wavelength-integrated emission in the wavelength interval 
(Ao + AAmin, A + AAmax) IS less than in the line's absence. Consequently, the P-Cygni 
absorption feature will be larger than the P-Cygni emission feature. (This ratio of feature 
size is not necessarily obtained if the line is not a pure resonance scattering line.) 

Now consider the toplighting case. First, imagine that there is no photosphere or line; 
there is just the radiating, optically transparent, spherical shell CSIR. The flux in any 
direction is 27csG'o at all wavelengths. Now add a hne to the atmosphere enclosed by the 
CSIR. The line has, in fact, no effect on the flux emission. The comoving- frame radiation 
field at the line center wavelength at any point inside the atmosphere is isotropic. The 
line scattering is isotropic. In the comoving frame, an isotropic field isotropically scattered 
is unchanged. Since the comoving frame radiation field is unchanged, the observer-frame 
radiation field is also unchanged. 

Now remove the line, but add a non-emitting (but, of course, opaque) photosphere. 
The flux in all directions at all wavelengths is now Ics (2Go — -^o)- The photosphere just 
acts as a net absorber. But if one now adds a line, the line will scatter some photons 
directed toward the photosphere into directions leading to escape. The result is that the 
wavelength-integrated flux in the wavelength interval (Aq + AAmin, A + AAmax) is increased 
by the addition of the line. 

Again one has to consider how the line scattering shifts the observer-frame wavelength. 
The non-emitting photosphere causes there to be a "cone of emptiness" in the radiation 
field converging at any point in the atmosphere. This reduces the line source function at 
every point in the atmosphere from the no-photosphere case. Thus, observer-directed beams 
in the limb region of the atmosphere that interact with the line must be reduced from the 
no-photosphere case. The closer the resonance point for a beam is to the z — location, 
the greater the effect of the "cone of emptiness" , and the weaker the emergent beam will 
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be. Beams that do not interact with the hne will be unchanged from the no-photosphere 
case. The upshot is that centered on Aq there will be flux deflcit relative to the continuum: 
an absorption feature around the line center wavelength. 

On the other hand, beams from the photodisk region (i.e., the non-limb part of the 
atmosphere on the near side of the photosphere) are enhanced by the line scattering. 
Without the line, the photodisk speciflc intensity is just /cs from the near hemisphere of the 
CSIR; the beams from the far hemisphere of the CSIR are occulted by the photosphere. But 
with the line there is scattering into the line-of-sight in the photodisk region. Consequently, 
there is an enhancement in flux over the continuum: i.e., a flux emission feature. Since 
the photodisk region is moving toward the observer this flux emission feature is blueshifted 
from Aq. 

One thus flnds that toplighting with a non-emitting photosphere gives rise to a flipped 
P-Cygni line with a blueshifted emission and an absorption centered near the line center 
wavelength. From the argument about wavelength-integrated flux above, we see that the 
blueshifted emission must be larger than the absorption trough centered near line center 
wavelength. 

In a general toplighting case the photosphere emits, and so there is a competition 
between P-Cygni and flipped P-Cygni line formation. Clearly in most supernovae the 
competition is won by P-Cygni line formation. For flipped P-Cygni line formation to win, 
E — Ics/I-ph must be greater than 1 (see § 4). 

6. Discussion 

The treatment of toplighting in supernova resonance-scattering line formation has 
been presented here in its simplest form for its heuristic value. For example, we have not 
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discussed how the tophghting might change the radial dependence of the hne optical depth 
r(r), and we have not considered the possibility that the CSIR reflects supernova hght 
back into the line-forming region. The toplighting version of the SYNOW code that we 
(Blaylock et al. 1999, in preparation) are using to analyze the spectra of SN 1998S allows 
for an angular dependence of the radiation from an optically thin circumstellar shell and 
for a wavelength dependence of the photospheric and CSIR specific intensities. Including 
toplighting allows us to obtain improved fits to the observed spectra of SN 1998S and shows 
that some of the P-Cygni hne profiles in the early-time SN 1998S may be fiipped. (It is not 
easy to be sure of the flipped proflles because of the numerous superimposed circumstellar 
features.) Lentz et al. (1999, in preparation) flnd that allowing for toplighting in detailed 
NLTE calculations also leads to improved flts. 

Among core-collapse events, SN 1998S has been uniquely well observed at ultraviolet 
wavelengths, but physically it is not an exceptional case. Line proflles in other circumstellar- 
interacting core-collapse events presumably also are affected by toplighting. But as we 
showed in § 4, only those supernovae whose monochromatic luminosities are strongly 
enhanced by circumstellar interaction will have line proflles that are strongly muted by 
toplighting. 

Examples of supernovae with relatively strong circumstellar interaction (as known 
fi:om relatively strong radio emission) are SN 1979C, SN 1980K and SN 1993J. AU these 
supernovae showed rather featureless UV spectra in the ~ 1800-2900 A region in comparison 
to supernovae known not to have had strong circumstellar interaction (Jeffery et al. 1994). 
Perhaps a UV-peaked tophghting continuum was muting the UV spectra of SN 1979C, 
SN 1980K and SN 1993J. 

Spectropolarimetry and nebular-phase hne proflles indicate that SN 1998S and its 
circumstellar shell were not spherically symmetric (Leonard et al. 1999) as has been 
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assumed here, and core-collapse events in general appear to be asymmetric (Wang et 
al. 1996). Eventually, asymmetric toplighting will have to be taken into account in spectrum 
calculations. 

This work was supported by NASA grant NAG5-3505, NASA grant GO-2563.001 to 
the SINS group from the Space Telescpe Science Institute, which is operated by AURA 
under NASA contract NAS 5-26555, and Middle Tennessee State University. 
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Fig. 1. — A simple model of supernova-line formation in the presence of toplighting. 



Fig. 2. — Supernova resonance-scattering line profiles in the presence of toplighting for 

-Res/ Rph = 2 
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